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Abstract 

Physical situations involving multiplicative noise arise generically in cosmology 
and field theory. In this paper, the focus is first on exact nonlinear Langevin equations, 
appropriate in a cosmological setting, for a system with one degree of freedom. The 
Langevin equations are derived using an appropriate time-dependent generalization 
of a model due to Zwanzig. These models are then extended to field theories and the 
generation of multiplicative noise in such a context is discussed. Important issues in 
both the cosmological and field theoretic cases are the fiuctuation-dissipation relations 
and the relaxation time scale. Of some importance in cosmology is the fact that 
multiplicative noise can substantially reduce the relaxation time. In the field theoretic 
context such a noise can lead to a significant enhancement in the nucleation rate of 
topological defects. 
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I. Introduction 

The study of noise in a field theoretic context has acquired ever increasing im- 
portance in the last few decades. More or less traditional applications such as phe- 
nomenological theories of spinodal decomposition have been augmented recently by 
several new research areas, e.g., stochastic quantization noise- induced transitions 
0, stochastic inflation and quantum decoherence [^]. The last two applications 
are of interest in cosmology. 

The objective here is to concentrate on the behavior of systems interacting with 
multiplicative noise in a few examples and to contrast it with what happens when the 
noise is additive. In the cosmological situation one also has to deal with the compli- 
cation that, since the Universe is expanding, there is a background time dependence 
of various parameters that would otherwise be fixed. 

Over the past two decades or so, the problem of understanding statistical field 
theory in a cosmological context has received much attention. One of the key issues 
here is to allow correctly for the time dependence of the background spacetime. Apart 
from any purely abstract motivation, this problem is of great practical interest as far 
as understanding the physics of the early Universe is concerned. Since not much is 
known about the early history of the Universe, approximations are often made which, 
although not justified rigorously, seem physically well motivated and even essential, 
if one wishes to obtain concrete results. 

A standard assumption is the idea that, at some level, the physical degrees of 
freedom of the Universe can be split into two coupled pieces, a "system" component, 
the evolution of which is to be followed in some detail, and a "bath" component, the 
detailed evolution of which is taken to be irrelevant. Examples of this picture include: 
(1) An understanding of inflation in terms of an infiaton field evolving in the rest of 
the Universe, which serves as an external environment [§-0], (2) the general notion 
of "coarse-graining" as a physical mechanism in terms of which to extract quantum 
decoherence, thus leading to a "quantum-to-classical" transition in the early Universe 
0, (3) the systematic development of statistical quantum field theory ||lOl, using a 
closed-time-path formalism to derive quantum dissipation and memory loss, (4) the 
intriguing, but not yet completely understood program of stochastic inflation 
originally proposed by Starobinsky [^. 

Much of the analysis has proceeded making use of essentially ad hoc Fokker-Planck 
equations, with an assumed bilinear interaction between the system and the bath. It 
is important to understand in what sense these models are physically reasonable, and, 
how relaxing the basic assumption of bilinearity changes the underlying physics, i.e., 
what happens when the noise in these models is multiplicative? 
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To understand some of these issues, it is instructive to consider models of systems 
interacting with baths consisting of time-dependent harmonic oscillators. In these 
models it is possible to derive exact, nonlocal Langevin equations |p. These mod- 
els are nonlinear, time-dependent generalizations of phenomenological models which 
have proven quite successful in other branches of physics, and, are well motivated 
phenomenologically, even if they are not derivable strictly from first principles. 

The analysis of Habib and Kandrup showed that (1) allowing for a nontrivial 
time-dependence necessarily induces qualitatively new effects like a mass (or fre- 
quency) "renormalization," even for the special case of bilinear couplings |]^; and, 
moreover, (2) allowing for nonlinearities in the system-environment coupling induces 
new effects aside from the usual "friction" term. For example, nonlinearities give rise 
to an additional renormalization of the system potential, and they lead to a "mem- 
ory" kernel which involves the state of the system. In the presence of nonlinearities, 
one's naive intuition is not completely lost: Even allowing for nonlinearities and time- 
dependent couplings in the interaction between the system and environment, for the 
case of time-independent oscillators (where the bath may still be viewed as being "at 
equilibrium" ) it is possible to derive a simple fluctuation-dissipation theorem [|12] . 
However, when the oscillators become time- dependent, the fluctuation-dissipation 
theorem is no longer exact. Nevertheless, the theorem is approximately valid as long 
as the coupling between the system and environment is dominated by modes of suffi- 
ciently short wavelength. In a cosmological setting this implies that, on scales short 
compared with the horizon length, it is still possible to speak of an approximate 
equilibrium and an approach towards that equilibrium. 

The examples treated here that are relevant in a field theoretic context refer to 
stochastic equations for such systems where the noise is a spacetime noise rather 
than the Gaussian, white noise encountered usually. By discretizing in space, it is 
possible in these cases to use oscillator models for the heat bath and to derive exact 
Langevin equations. Since these equations are at the heart of one approach to finite 
temperature simulations of field theories, it is useful to have a simple (indeed trivial) 
derivation for them. The two specific cases considered here will be a double-well 
0^ theory in 1 + 1-dimensions and the nonlinear sigma model in 1 + 1-dimensions. 
The last example is an interesting one: even though the coupling to the heat bath 
is bilinear, the noise is multiplicative. This happens because the stochastic forcing 
term must be consistent with a constraint. 

Section II will be devoted to the oscillator model and its generalizations. The equa- 
tions of motion will be derived and the conditions under which a Langevin description 
arises will be discussed. Particular attention will focus on fiuctuation-dissipation re- 
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lations, time dependent baths, multiple degree of freedom systems, and nonlinear 
system-bath couplings. Section III will deal with the modification of the relaxation 
time scale by multiplicative noise. The method of Lindenberg and Seshadri [|T3] (first 
suggested by Stratonovich [|^) will be applied to a cosmo logical model problem. It 
will be shown that substantial changes in the relaxation rate can occur. Two examples 
studied recently provide numerical evidence: a harmonic oscillator driven to equilib- 
rium by quadratic multiplicative noise and the nucleation of kinks in the double-well 
0^ theory in 1 + 1-dimensions. In the first case, the enhancement of the relaxation 
rate predicted by the theory is consistent with the simulations, while in the second 
case an enhancement of the kink nucleation rate is observed (this process has not 
yet been investigated analytically). Finally, Section IV will offer some conclusions, 
unanswered questions, and speculations regarding future results. 
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II. Multiplicative Noise, Time Dependence, and the Oscillator Heat Bath 

The material at the beginning of this section borrows heavily from Ref. ^ where 
the reader will find a rather more detailed discussion of the derivation presented here. 

The idea behind the oscillator heat bath model is relatively simple: given some 
arbitrary composite Hamiltonian one splits it into three pieces: the system, the bath, 
and the interaction. If this formal split is to be useful, the system must be "small" in 
some sense when compared to the bath or environment. As a practical matter what 
this means is that, as far as the system is concerned, the full Hamiltonian 

Ht = Hs + Hb + Hi (1) 

is well approximated by 

H = Hs + 6Hb + SHi (2) 
where SHb is the Hamiltonian for a collection of harmonic oscillators and SHj an 



interaction Hamiltonian linear in the oscillator variables qa [T^ ||T6|. Basically, one 
assumes that, in the absence of any coupling with the system, the environment is 
characterized by some fixed, possibly time-dependent, solution. Allowing for a weak 
coupling with the system (weak in the sense that each bath mode is only changed 
marginally), one then identifies the g^i's as perturbed variables, i.e., degrees of freedom 
defined relative to the background solution (e.g., phonons). This has two implications: 
(1) the environment can be visualized as a collection of oscillators with (possibly time- 
dependent) frequencies, so that 6Hb is quadratic in bath variables qa, and, (2) because 
the interaction of the environment with the system is assumed to be weak, 6Hj must 
be linear in the g^'s. Note that the system is not necessarily weakly altered, so the 
interaction SHj does not have to be linear in the system variable x. In principle one 
can proceed without imposing any restrictions on the form of the system Hamiltonian 
Hs. 

With the above set of assumptions, one can write 

Hs = ^V^ + Vren{x,t), (3) 



SHb , 



^E[pa+^aW?a]. (4) 

^ A 

and, in terms of relatively arbitrary functions Ta, 



SHi = -Y.Q^{t)rA{x,t)qA. (5) 

A 

It is, however, convenient to rewrite H in the manifestly positive form 

H = ^v' + uix, + ^ E {pI + (i) [qA - r^(x, t)f} , (6) 
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where in terms of the "renormahzed" potential K 



reni 



U{x,t) = Vren - ^ J] fi^ (t)!^ (o;, t) . (7) 
^ A 

Couphngs of a system to some environment can induce finite and stochastic renor- 



mahzations in the system potential, although this is not always so |T6[ . In this paper, 
the words "system potential" will always refer explicitly to the renormahzed system 
potential. 



Finally, as emphasized, e.g., by Caldeira and Leggett |T^, it should be stressed 
that this is more than simply a toy model. This form of the Hamiltonian generally 
provides a correct description for any system which is only weakly coupled to its 
surroundings (Although there are caveats relating to the spectral distribution of the 
modes: one has to be sure that there are no negative frequency, unstable modes. 
This is not a trivial issue when the full Hamiltonian is complicated and not well 
understood). To facilitate a concrete calculation, this Hamiltonian need only be 
supplemented by two inputs, namely the spectral distribution of the environmental 
modes and the form of the coupling to the system. For a general physical problem, 
these may either be extracted from experimental data or derived from theoretical 
considerations (of course, this may not turn out to be simple in practice). 

The equations of motion generated from the Hamiltonian (^ clearly take the forms 



X 
V 



-— U{x, t) + Y: W [U - TAix, t)] Q^r^ix, t), 



Qa = Pa, 

Pa = -n^[qA-TA{x,t)], (8) 

where an overdot denotes a time derivative d/dt. Following Ref. one can show 
that these equations collapse to the single nonlocal equation 

V = - / ds [Kit, s)v{s) + M(t, s)x{s)] + F,(t), (9) 

ox Jo 

where, in terms of the quantities 

AA{s,t) = (y.iT\s)x-\s)^ (E7?^Wx"-^(t)) (10) 

and 

BA{s,t) = -—AA{t,s), (11) 
n OS 

the "memory kernels" K{t, s) and M{t, s) are 

K{t,s) = Y.^Iit)Ms,t)WA{s,t), (12) 

A 

M{t,s) = J2^Iit)BA{s,t)WA{s,t). (13) 

A 
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Here one has assumed a polynomial system-bath coupling 

N 1 

TA{x,t) = Y.-l^A\t)x^ (14) 

and the force Fs is 

Fsit) = J2n^{t)J27^r\t)x"^-' {MO) -TAix,o)]CAit) +pAio)SAit)} , (is) 

A rn 

where 5'^(t) and Cyi(t) denote two linearly independent solutions to the homogeneous 
oscillator equation 

EA + n^mA = 0, (16) 

satisfying the initial conditions Ca(0) = Sa{0) = 1 and Sa{0) = Ca{0) = 0. 

Equation (|^) does not at first sight appear to be an ordinary Langevin equa- 
tion. As a first step to see how this may be so, it is possible to prove a general- 
ized fluctuation-dissipation theorem, even if Ta is a nonlinear function of x and/or 
explicitly time-dependent |P](as long as the oscillator frequencies Qa are not time- 
dependent). 

Consider an ensemble of initial conditions for which the first moments vanish 
identically, i.e., 

(Qa(0)) = (pa(0)) = 0, (17) 

with Qa{0) = qA{0)—TA{x{0), 0), and where the second moments are initially thermal, 
so that 

(pa(o)pb(o)) = ^]A(o)^^B(o) (gA(o)gB(o)) = ksTSAB, (is) 

where the angular brackets denote an initial ensemble average. Then 

(Fsit)) = (19) 

and 

{Fs{h)F,{t2)) = kBTK{t,,t2), (20) 

thereby identifying Fs(t) as a noise and proving a generalized fluctuation-dissipation 
theorem linking the noise autocorrelator with the "viscosity kernel" K{t,s). Equa- 
tion (^ can now be viewed as a nonlinear, nonlocal Langevin equation. It is im- 
portant to note that for arbitrary time dependence of the oscillator frequencies, the 
fluctuation-dissipation theorem does not hold. In the cosmological context, an ap- 
proximate fluctuation-dissipation relation is valid as long as the system time scale is 
small compared to the expansion time scale of the Universe. 

The Langevin equation (|^) reduces to a well known form in an appropriate limit. 
If one neglects nonlinearities in the coupling between system and bath, assuming that 
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oc X, one immediately recovers a special model considered in Ref. ||^. With the 
further assumption that Qa and Ta are independent of time, one is reduced to the 
well known independent oscillator (10) model |T^. It is thus possible to address 
systematically the question of how the incorporation of nonlinearities and/or time- 
dependences leads to systematic changes in the Langevin equation derived for that 
original model. Ref. ||^ discusses such issues as time scale separations and Markov 
limits (when the nonlocal in time equation becomes local). Here it will be simply 
assumed that it is indeed possible to write (P) as a local, albeit nonlinear Langevin 
equation. Once this can be done, it is immediate to write the Fokker-Planck equation 
(for the time evolution of the phase space distribution function) that corresponds to 
a particular Langevin equation. 

Can this simple approach be extended to continuous systems, e.g., field theories? 
When discretized, classical field theories can be thought of as interacting continuous 
spin systems. In this case, the model works more or less as usual except that the 
system has a very large number of degrees of freedom and the effective noise turns 
out to be a "spacetime noise" rather than the one considered so far. The Langevin 
equations derived in this way are useful for thermal simulations as an alternative 
approach to Monte Carlo techniques. In what follows no explicit time dependences 
will be assumed. 

As a relevant example of a field theoretic problem, a double-well X(f)^ theory in 
1 + 1-dimensions is particularly convenient. This theory admits stationary wave 
solutions ("kinks") and the statistical mechanics of these objects has recently been 
carefully investigated via Langevin simulations []18 |. 

The Lagrangian for this theory is 

L = l(a,$)2 - i(9.$)2 + - 1a$^ (21) 

In order to derive the appropriate Langevin equation for this theory it is convenient 
to discretize in space and use the Hamiltonian formulation. The discrete Hamiltonian 



Pi - m $i + -X<l>i + - + 



(22) 



2*2V A / 2V A 
gives rise to the equations of motion 

92$, = A^$, + $,(1 - $,') (23) 



where A is the lattice spacing and the spatial lattice Laplacian 



It is clear that H[) can be interpreted as describing a continuous spin system with spins 
$j (at the lattice site i) interacting via nearest neighbor interactions. To drive this 
system to equilibrium the idea is now to couple each spin to an independent oscillator 
heat bath: there will now be as many "baths" as there are spins, all with the same 
initial conditions in the sense that the averages (|l^ and (|l^) hold for each collection 
of oscillators. By following the derivation given above, but for each individual spin 
this time, it is easy to derive the Langevin equation (written, for simplicity, for a 
linear coupling to the heat bath) 

92$, = A^$, - r]dt<^, + $,(1 - $,2) + F, (25) 

where the noise auto-correlation 



{F,{t)F,{s)) = 2'nkBT5{t - s)(5,,/A). (26) 

Note that by taking independent heat baths for the spins we have enforced the van- 
ishing of the noise cross- correlation. Since the fluctuation-dissipation relation is still 
valid, we know that that the system of Langevin equations (^) will drive the system 



to thermal equilibrium. The continuum limit of (|25|) is the Langevin equation 

dl^ = - Vdt<^ + $(1 - $') + F{x, t) (27) 
and in this limit the noise correlator 



(F(x, t)F{y, s)) = 2r,kBT5{t - s)5{x - y). (28) 

describes a white, "spacetime" noise. 

One way of introducing multiplicative noise is through the imposition of con- 
straints on the system variables. One enforces the noise to act in such a way that 
the reduced dynamics for the system respects the constraint imposed via a Lagrange 
multiplier. This process turns what would have been additive noise into multiplica- 
tive noise: even if the coupling to the environment is linear in the Hamiltonian, the 
effective Langevin equation contains multiplicative noise. 

As a representative example, consider the 0(3) nonlinear sigma model in 1 + 1- 
dimensions. The full Hamiltonian is 

H = J dx ^{p-p + d^h ■ d^h) , (29) 

with the constraint, 

h-n = l. (30) 
9 



To implement the constraint, one introduces a Lagrange multiplier A and writes the 
Hamiltonian as 

H = J dx -{p ■ p + dxh ■ dxh + \{h ■ n — 1)) . (31) 

The equation of motion is then 

(-4 + 5L)^ + An = (32) 

The Lagrange multiplier can be eliminated from this equation by taking the dot 
product of it with h. This procedure yields 

{-dl + dlx)n - {n ■ {-dl + aL)n)n = 0. (33) 

One can now discretize the Hamiltonian in much the same way as the (^"^ theory 
considered above. Following that analysis one linearly couples oscillator heat baths to 
each spin component at each lattice point and derives in the usual way, the Langevin 
equation 

din - dl^h - An + rjdtfi + F = (34) 
Again, by eliminating the Lagrange multiplier, one obtains 

din — dl^h + (n ■ dl^n)n + {pth ■ dtn)n + rj^dth — {n ■ dth)n) + F — {n- F)h = 0, (35) 

a nonlinear Langevin equation with multiplicative noise. This equation (supple- 
mented with a symmetry breaking term) has been used in a numerical study of 



thermal sphalerons in the nonlinear sigma model [jT9| where it has been verified that 
it takes the system to the correct equilibrium distribution while maintaining the con- 
straint to a reasonable level of accuracy (1 part in 10^, the error due to time-stepping 
inaccuracies). This way of driving the sigma model to equilibrium is better than 
using variables which explicitly maintain the constraint (for example, the polar an- 
gles) since there are no coordinate singularities to cause dynamic range problems in 
numerical simulations. 



10 



III. Relaxation with Multiplicative Noise 

In the previous section multiplicative noise entered in the 10 model in two ways, 
first through a nonlinear coupling to the environment and second, through the impo- 
sition of a constraint on the system variables (which effectively introduces a nonlinear 
coupling). Systems subjected to noise of this kind can have striking departures in 
their dynamical behavior when compared to the case of additive noise. Here the fo- 
cus will be on the process of relaxation to equilibrium: Even though multiplicative 
noise takes the system to the usual thermal equilibrium state (when the fluctuation- 
dissipation theorem is valid), the relaxation time scale can be quite different from 
that for additive noise. In particular, it can be much smaller. 

As an illustrative example, consider the cosmologically relevant 10 model of a 
"system" variable interacting with scalar "radiation" in a spatially fiat, Friedmann- 
Robert son- Walker universe (discussed more fully in Ref. [§, from where this example 
is taken) 

H = \v^ + U{x, + W + ^liv) [Qa - r^(x, v)]}' , (36) 
where, in terms of constants A, ^, and o", 

[/(x, r/) = \Kri)x' + U{ri)r' + \a{j])x\ (37) 

and the coupling 

YA{x,ri)=-i'^{r])x^-if{r])x^. (38) 

The frequencies VLa satisfy 

^l{ji)=ujl-{\-^i)- (39) 



where a prime ' denotes a conformal time derivative djdri and a is the scale factor. 
The conformal time ri is related to the cosmic time t via dri = a'^dt. 

If the scale factor a has a simple power law time dependence a = t^, with p > 1/2, 
it follows that 

aoc{7]~r]oy/^^-P\ (40) 
where tjq denotes an integration constant, so that 



where 



"^-(i-eof^ (42) 
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is intrinsically positive when p > 1/2 and 1 — 6^ > 0. The quantities denote 
eigenvalues of the spatial Laplacian. 

To proceed further, one obviously needs to determine the mode functions for 
the time-dependent frequencies. The equation 

,2 



iv-Vo) . 



is solved by a general 



(43) 



(44) 



where denotes an arbitrary solution to Bessel's equation of order 



= (1 - GO 



p{2p-l) 1 



2 1 

4 4' 



(45) 



(l-p)2 

For the special cases ^ = 1/6 (conformal coupling) and/or p — 1/2 (a Universe 
dominated by electromagnetic radiation), /i — 1/2 and the solutions reduce to 
sines and cosines. For these particular values, the bath Hamiltonian 5Hb is time- 
independent in the conformal frame. This implies that the physical frequencies are 
simply red-shifted uniformly as the Universe expands. 

Even when the special conditions given above do not hold, as long as the system 
time scale is small compared to the expansion time scale, one recovers a relatively 
simple nonlocal Langevin equation, for which the only explicit time dependences are 
in the potential U and the couphngs between the system and the environment. It 
thus follows that, in this approximation, a fluctuation-dissipation theorem holds, so 
that one would anticipate an evolution towards some steady state solution at late 
times. The fluctuation-dissipation theorem is of course exact when ^ = 1/6 and/or 
p=l/2. 

If a Markov approximation is valid, the nonlocal Langevin equation can be well 
approximated by a local equation of the form 



dU 

v' = - M{v)x{ri) - IC{v)v{v) + Fs{v): 



(46) 



where, for an initial thermal ensemble. 



(F,(77)F,(s)) = 2kBT)C{x,r])6Div - s) 



(47) 



This local Langevin equation leads immediately to a Fokker-Planck equation of the 
form 



df] dx dv 



dU 
dx 



Mx - JCvj f 
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ksTK^^Q. (48) 



Translated back into the physical frame, the Langevin equation (56) becomes 



V 



1 dU 



m{t)x 



a 



3- + K(t) 

a 



a^dX 

where the overdot denotes differentiation with respect to cosmic time, and 

X 

a 



X -- 
a 

M 



X = 
V = 

m = 

Mit) = 

Fs(t) = 
The noise autocorrelator is now 

(F.(t)F.(t')) 



V 



a 



Fs 
a3 



+ K(t) 



+ 



+ 



(49) 

(50) 
(51) 
(52) 
(53) 
(54) 

(55) 



Note that, in the physical frame, there are two sources of damping, namely the 
viscosity oc K.V and the cosmological frame- dragging oc H = a/a. The corresponding 
Fokker-Planck equation is now 



d d 
dt ' SJ^^^^^+^y 



1 dU 
'^dX 



M{t)X- {3- + K{t))V 



f l-^K(t)^ 



0. 



(56) 

Note that the conformal temperature T is related to the physical temperature Tph by 



ph 



T/a. 



The analysis of the approach to equilibrium now follows that of Lindenberg and 
Seshadri [|1^] based on earlier work by Stratonovich [|14|. The basic idea is that for 
sufficiently weak damping the time scale on which the system coordinate changes is 
much faster than the time scale on which the energy evolves. This allows one to 
average over the fast coordinates and convert the original Langevin equation, written 
in its Fokker-Planck representation, to a one-variable Fokker-Planck equation for the 
energy distribution function (hence the name, "energy-envelope method"). 

As a first step one implements a change of variables from (x, v) to (x, E), where, 
E = +U{x,ri), to obtain a new Fokker-Planck equation for the distribution 
W{x,E,t) satisfying 

/(x, V, rj)dxdv = W{x, E, rj)dxdE. (57) 

To the extent that the energy is approximately conserved during a single oscillation 
of the system, one can assume that 

W{x, E, r]) = {2$'(E) [E - U{x, r])]^^'^}-^Wi{E, r]) 
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(58) 



where 

= J dx[E- U{x, r])]^^^ (59) 

and a prime now denotes a d/dE derivative. Here the integration extends over the 
values of x along the unperturbed orbit associated with E. Note that the prefactor 
of Wi is simply the relative amount of time that, for fixed energy, the system spends 
at each point x. 

By integrating the Fokker-Planck equation for W{x, E, r]) over x, one obtains the 
desired equation for Wi{E,ri), which takes the form 

d 



dr] 



W,{E,r]) 



'd_ 

dE 



1 



X'{E) 

+\2{^{E)-kBT^\E))}] 



{Ao {^E) - kBT^'{E)) + 2Ai {x{E) - kj.Tx'iE)) 



+kBT 



dE^ 



{Xo^E) + 2Xix{E) + \2^iE)} 



where 



and 



X'{E) 

X{E) = J dxx [E-U{x,r])] 
^(E) = J dxx^ [E-U{X,T])] 



H^i(E,r/), 



1/2 



1/2 



(60) 



(61) 



(62) 



Unfortunately, for a generic potential U the functions $(-£'), x{E), and "^{E) can- 
not be evaluated analytically, so that one cannot realize the right hand side of ( |62D 
explicitly in terms of simple functions of E. For a quartic potential, for example, 
these functions can only be expressed as elliptic integrals, which must be evaluated 
numerically. In the limit when the energy E is small and the system is oscillating 
about a local minimum of U one can proceed analytically by evaluating the orbit 
integrals, assuming that the system is effectively evolving in a simple harmonic oscil- 
lator potential. By doing so, it is straightforward to derive from the Fokker-Planck 
equation a transport equation involving the time derivative of the first energy moment 



{E{r])) = j dE E Wi{E,7]). 



(63) 



First suppose that S is positive. Then the moment equation takes the form 
I- (E(,)) = kBTX, - f Ao - ^1 (i?(,)) - ^ lE\r,) 



dr] 



(64) 



^0 / ^0 

Note that, because of the reflection symmetry x —>■ —x for the potential U, the 
functions x = x' = 0, so that the contributions involving Ai vanish identically. 
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Unfortunately, this equation still cannot be solved exactly for {E{ri)), as it involves 
the unknown function {E'^{ri)). To obtain a formula for {£"^{7])), one must consider 
the second moment equation, which in turn relates d {E'^{ri)) /dr] to the third moment 
{E^{ri)). In the spirit of (say) the BBGKY hierarchy, one requires a truncation 
approximation. 

As in Refs. []T^ and ||TB[, suppose that 



(E'it)) ^ K{E{t))\ (65) 

with K = 2. One knows that, when the system is at equilibrium, with energy E = ksT, 
this equation is satisfied identically for k = 2, and one might expect on physical 
grounds that, before the system is "at equilibrium," the energy distribution will be 



narrower and k < 2. As emphasized by Lindenberg and Seshadri this truncation 
approximation thus leads to an upper limit on the time scale on which the system 
"equilibrates" with the bath. Given this truncation, one can immediately write down 
the solution 



(P( ^\ _ ''bT {Eq + AksT) - AkT {ksT - Ep) exp {- [{A + 1)/A] Xpf]} 

^ {Eo + AkBT)-{kBT-E,)e^^{-[{A+l)/A]\,r^] ' ^ ^ 

where 

A ^ (67) 

In the limit that A2 ^ 0, the system approaches an "equilibrium" with {E) = ksT 
on a time scale Ir ~ Aq ^. If A2 7^ 0, the system still evolves towards an equilibrium 
with {E) = ksT, but the time scale can be quite different. In fact, in the limit 
that Ao — > 0, 

tR^LU^yiX^ksT). (68) 

Therefore, when the nonlinear coupling is sufficiently strong, the system may be 
driven towards equilibrium, not by the ordinary additive noise associated with the 
linear coupling, but primarily by the multiplicative noise associated with the nonlinear 
coupling. 

The only point that remains to be checked is that one is still assuming, as is 
implicit in this "envelope" approximation, that the time scale ujq^ is much shorter 
than the damping time. This, however, is clearly the case when Aq and A2 are not 
too large. Indeed, one verifies that (a) the weak damping approximation is legitimate 
but (b) multiplicative noise dominates the evolution towards an equilibrium whenever 

m 

^ « ^ « 1. (69) 
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Consider now the case S < and the system oscillating about one of the two min- 
ima of the potential xo^O. (This is the case relevant to first order phase transitions.) 
Here the terms involving x and its energy derivative do not vanish, and the formulae 
for $ and \1/ acquire additional terms involving the location Xq of the new minimum. 
However, one still recovers a relatively simple exact equation for d (E) /drj. In this 
case, the moment equation is 

A iEi„)) ^ tsTL - (i - ^) - ^ {S^ir,)) . (70) 

where 

S = E-U{xo) (71) 
denotes the system energy defined relative to the minimum of the potential, and 

L = Ao + XiXq + Aoa;Q^ (72) 

plays the role of a "dressed" Aq. Note in particular that, if |xo| is large, as will be 
the case when S <C A, the dressed L can be much larger than Aq. This implies that, 
in this case, the nonlinear couplings can reduce the overall equilibration time both 
through the introduction of a new term oc A2 (-E^) and through an increase in the 
effective linear coupling. 

In any case, as long as U{xo,t]) is only slowly varying in time, the derivative 
d {Eirf)) /drj can be replaced by d {£{ri)) /drj. And, to the extent that the coefficients 
L and A2 may be approximated as time-independent, (E9|) can be solved analytically. 



The result is an expression identical to (^21) , except that E is replaced by the shifted 
£ = E — Uq and 

Numerical tests of these results are not difficult to carry out. A simple test is to 
take a harmonic oscillator and couple it quadratically to noise. Simulations show, that 
in the weak coupling regime, the relaxation to equilibrium is indeed much faster than 
for the case of linear coupling and that the relaxation time scale is in agreement with 
(|68|) pO| . In the case of strong coupling to the heat bath, the energy-envelope method 
is no longer applicable. Numerical simulations indicate that the relaxation time scales 
for both the additive and multiplicative noises are similar at large coupling. 

An interesting example that corresponds to the case S < occurs in the problem of 
determining the nucleation rate of kinks in a double- well 0^ theory in 1 + 1- dimensions. 
One begins with the system in one vacuum and counts kinks as they nucleate (as pieces 
of the field configuration "hop" to the other side). At low temperatures compared to 
the kink energy, the nucleation rate is very low (and is a significant problem when 
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carrying out numerical simulations). However, in the presence of multiplicative noise 
the rate can be significantly enhanced by following an argument exactly analogous 
to the one given here for the particle system. While the theory has not yet been 
completely worked out, numerical simulations do show a significant increase in the 



nucleation rate in the case of nonlinear noise 20 
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IV. Conclusion 

The objective here has been to motivate and derive exact Langevin equations for 
particle systems and field theories using an oscillator model for the heat bath. Multi- 
plicative noise is seen to arise in these models as a consequence of nonlinear coupling 
to the heat bath, a situation particularly natural when dealing with systems with 
constraints. A particularly interesting consequence of nonlinear noise is the drastic 
change possible in relaxation time scales: a simple analytic estimate is possible in 
the case of weak damping using the energy-envelope method. Numerical simulations 
have been carried out to confirm the theoretical prediction. While not yet complete, 
they are in agreement with the theory so far. 

There are several interesting applications of this approach. Langevin simulations 
for determining the rate of sphaleron transitions at finite temperature in the nonlinear 
sigma model have recently been carried out using the nonlinear Langevin equation 
derived here. Such simulations are also easy to carry out for spin chains using es- 
sentially the same technology. The problem of the dependence of the nucleation rate 
of kinks on the nature of the noise (additive or multiplicative) is another interesting 
issue. There is clear numerical evidence that the nucleation rate is speeded up by 
multiplicative noise but a complete analytic theory has still to be worked out. Work 
in this direction is in progress. 

One of the original motivations for this work was to study the onset of inflation in 
such a way as to answer the question first raised by Mazenko, Unruh, and Wald ||21|| : 
Does a homogeneous slow-roll really take place or does domain formation happen 
so quickly that infiation is ruled out? This question was tackled by Albrecht and 



Brandenberger |]22[ who argued that for a sufficiently small value of the self-coupling 
for the infiaton field the scenario remains viable. Recently, the effect of fiuctuations 
was studied in Ref. via an essentially ad hoc approach, the results of which did 
not weaken the estimates of Ref. |Q. Now, when one deals with an infiaton field 
with a (small) quartic self-coupling, the effective equations for the mean field contain 
a (weak) nonlinear damping term. It is then reasonable to ask what happens when 
the relaxation to the true minimum is controlled by these terms. The work presented 
here suggests that these terms can dramatically change this time scale. Whether the 
change is of practical relevance can only be determined by further calculations, both 
analytic and numerical. 
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